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The measured partial widths of the semileptonic decays D — > K£v and D — > niv can be combined 
with the form factors calculated on the lattice to extract the CKM matrix elements |V CJ | and \V cc i\. 
The lattice calculations can be checked by comparing the form factor shapes from the lattice and 
experiment. We have generated a sizable data set by using heavy clover quarks with the Fermilab 
interpretation for charm and asqtad staggered light quarks on 2+1 flavor MILC ensembles with 
lattice spacings of approximately 0.12, 0.09, 0.06, and 0.045 fm. Preliminary fits to staggered 
chiral perturbation theory suggest that we can reduce the uncertainties in the form factors at 
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1. Introduction 

Performing a global fit under the assumption of a unitary CKM matrix ^ yields precise values 
for \V CS \ and \V c d\ M{. If new physics in flavor introduces deviations from unitarity, however, values 
of the CKM matrix elements from direct determinations will in general differ from those predicted 
by the global fit. Moreover, improving tests of second row and column unitarity requires improved 
direct determinations of \V CS \ 

In the limit of massless leptons, the rates for the semileptonic decays D — > K(n)lv become 1 



dT(D ^ K(Tl)lv) = G 2 F u |3|v A 2 U D^K(n) ( ^2 
dq 2 24n 



where q 2 = (p D — Pk(tz)) 2 ^ Pr(w) * s tne momentum of the K(ri) in the rest frame of the D, and 
j.d-^k(h) j s defined in terms of the hadronic matrix element of the current = isy^c (idy^c): 

(2 2 \ 2 2 
m D~ m K(n) \ , j.D^K{K) t 2 \ mD ~ m K{x) 
Pd + Pk(k) T^y] + fo (*) 2 — Qv- 
q J, q 

(1.2) 

Given the normalization of the form factors y^^W f rom lattice QCD, the CKM matrix ele- 
ments |V C4 .(^) | can be extracted from experimental measurements of the branching fractions. 

Agreement with the Standard Model values provides important validation of our methods, 
which we also use to calculate the form factors for B — > %lv [Q] and B — > K£ + £~ Both decays 
are central in searches for new physics; the former allows direct extraction of \V u b\, while the 
latter is loop-suppressed in the Standard Model. New physics seems unlikely to enter the tree-level 
decays D — > K(n)lv before the B decays. The Fermilab method applies to charm and bottom, so 
consistency between our D form factors and the normalizations implied by the global fit is direct 
evidence of our ability to precisely extract the B form factors. 

2. Method 

For the up, down, and strange quarks we use the asqtad-improved staggered action [^], for the 
charm quark we use the clover action with the Fermilab interpretation [RJ], and for the gluons we 
use a one-loop Symanzik improved gauge action [Js[] . We set the scale with f„, tune the masses of 
the light quarks using the experimental values of m % and thk , and tune the hopping parameter of 
the charm quark with the spin-averaged kinetic mass of the D s . 

Table [j] summarizes our data set. We vary the valence light-quark masses on each ensemble 
from near the tuned strange mass m s down to ~ 0. \m s , and the lattice spacings from w 0. 12 fm to 
w 0.045 fm. To increase statistics and reduce autocorrelations, we average over four source times 
and randomize the source spatial locations. 

For calculations in the rest frame of the D and in heavy-meson chiral perturbation theory, the 
hadronic matrix elements are conveniently parametrized by form factors f± and /j| : 



(K(K) \V, |D) = V2m~ D fv^W iq 2 } + p ^f^) (^) 



(2.1) 



For decays to the % , there is an additional isospin factor of 1 /2 on the right-hand side of Eq. (1.1) 
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~ a (fm) 


I? x n t 


x * corn 


^src 




am/ / flWZ iS - 




' c 


0.12 


20 3 x64 


2052 
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0.02/0.05 


{0.005, 0.007, 


0.1259 




20 3 x64 


2259 
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4 


0.01/0.05 


0.01, 0.02, 


0.1254 




20 3 x64 


2110 


4 
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0.007/0.05 


0.03, 0.0415, 


0.1254 




24 3 x 64 


2099 
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0.005/0.05 


0.05, 0.0349} 


0.1254 


0.09 


28 3 x96 


1996 
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0.0124/0.031 


{0.0031, 0.0047, 


0.1277 




28 3 x96 


1931 


4 
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0.0062/0.031 


0.0062, 


0.1276 




32 3 x96 


984 


4 
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0.00465/0.031 


0.0093, 0.0124, 


0.1275 




40 3 x 96 


1015 
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4 


0.0031/0.031 


0.031, 0.0261} 


0.1275 




64 3 x96 


791 


4 


4 


0.00155/0.031 




0.1275 


0.06 


48 3 xl44 


593 


4 


4 


0.0072/0.018 


{0.0018, 0.0025, 


0.1295 




48 3 x 144 


673 


8 


4 


0.0036/0.018 


0.0036, 


0.1296 




56 3 xl44 


801 


4 
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0.0025/0.018 


0.0054, 0.0072, 


0.1296 




64 3 x 144 


827 


4 


4 


0.0018/0.018 


0.0160, 0.0188} 


0.1296 


0.045 


64 3 x 192 


801 


4 


4 


0.0028/0.014 


{0.0018, 0.0028, 


0.1310 



0.0040, 0.0056, 0.0084, 0.0160, 0.0130} 



Table 1: Data on the 2+1 flavor asqtad staggered MILC ensembles for various valence masses and source 
times. The columns are, respectively, the lattice spacing, lattice dimensions, number of configurations, 
number of source times, number of 3-point sink times, light/strange sea-quark masses, valence light-quark 
masses, and charm hopping parameter. The analysis to date includes the full QCD points; we are considering 
generating data at additional source times. 



where v = po/mD, and p± = Pk(ti) ~ {Pk{ti) ' v ) v - fx and f\\ can be extracted from correlator ratios. 
We consider 

igP-*(«)/ fT .J-\- 1 C^ K{K \t,T;p K{7t) ) 2E K[K) 

where tj> K{7l)fl = (1, p K (n)), E K(n) = ( m D + m \(n) ~ q 2 )/{^m D ), andC 3 , C 2 are averages of corre- 
lators constructed to eliminate oscillations from opposite-parity states [Q]. T and t are respectively 
the source-sink separation and current insertion time in the vector-current 3-points. 

We use local operators for the K(n) 2-points, smear the D interpolators with a charmonium 
wavef unction, and construct the currents out of light staggered and heavy clover fields [§]. For 
insertion times t far from 3-point source and sink, the ratios plateau to the form factors: 

R 3 Q ~ f« y ' and R 3 i ~ f x w for 1 < t <C T. (2.3) 

The averages C3 require raw 3-points at successive source-sink separations T and T + 1 ; to min- 
imize our errors as a function of momentum p^), we generate the 3-point correlators at two 
physical separations on each ensemble and for each set of valence quark masses JTOj]. 

We inject the 3-points and K(n) 2-points with momenta ap K ^/(27i/L) = (0,0,0), (1,0,0), 
(1,1,0), (1,1,1), and (2,0,0) (and permutations and negatives of these components). We aver- 
age the correlators over equivalent momenta (up to axis interchange) and have checked that the 
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wavefunction overlap factors of the K(ti) 2-points are independent of momentum. We then replace 
C2^\pk(k)) with the less noisy C^^(O) in the ratios of Eq. (22) and use 




C2 W (t;0)Ci(T 

We have checked that our data obeys the continuum dispersion relation, and we substitute this 
relation for E K i % \ in the above ratio. 

To extract the masses of the K{%) and D from the 2-point correlators, we fit to sums of ex- 
ponentials with oscillating terms to account for contributions from opposite-parity states. We use 
the masses to construct the ratios, and we propagate the errors via 500 bootstraps. To extract the 
plateaus we fit the ratios; varying the fit function and time intervals does not significantly change 
the results. 

The vector currents undergo renormalization. We match to the continuum by writing the renor- 
malization factors as products of the degenerate-mass vector-current normalization factors, which 
we compute nonperturbatively, and correction factors whose deviations from one are perturbatively 



calculable [11 j: 



(K(K)\V,\D) =Z^(K(x)\V l *\D), Z c v ^=p c v f^Z^ dd \ (2.5) 

We blind the analysis by introducing an offset in p. 

After operator renormalization, the results for the form factors on all ensembles and for all 
combinations of valence masses and momenta are simultaneously fit to staggered chiral perturba- 
tion theory (S^PT) Q. Fits to SU(3) S^PT are shown in Fig. [I[ The data are for the D -> % decay 



on a subset of the ensembles of Table [[]: the four coarse (a « 0.12 fm) ensembles, four of the fine 
(a ?» 0.09 fm) ensembles (excluding the 0.05ra s ensemble), and three of the superfine (a « 0.06 
fm) ensembles (excluding the 0.14m s ensemble). Data with momenta up through = (1, 1,0) are 
included in the fits. 2 The fit function includes leading chiral logarithms (NLO loops) and analytic 
terms through NNLO. 

The SU(3) S^PT fits for jf^* and ff~* K are good, much better than SU(3) S^PT fits for 
jf-+* and f^ K . As exemplified in Fig. | fits of data to SU(3) S^PT are marginal in 

quality. We understand this behavior in terms of the features of the #PT description of the energy 
dependence of the form factors. 

The pole from the resonance dominates the energy dependence of f±, but the ^PT expressions 
for /|| contain no pole, and the fits suffer. The fit quality for fu may also reflect the limitations 
of #PT. For energies comparable to the chiral symmetry breaking scale, E K ^\ ~ A^sb. we expect 
the #PT description of the energy dependence to break down. The largest-momentum data in 
fits to date have p^-(^) = (1,1,0), which corresponds to E K ~ A^sb; on the coarse ensembles, 
Pn: = (1,1,1) corresponds to energies E n > A^sb- 

SU(3) S#PT fits for fn*~* K support the idea that the #PT description of the energy dependence 
of /|| is breaking down. For the same momenta, the energies Ex are larger than the energies E n , the 



2 Hereafter all momenta are given in units of 2n/(aL) unless otherwise specified. 
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#PT expansion breaks down for smaller momenta, and fits for fP should be worse than fits for 
fP-+ x . In fact, SU(3) S^PT fits for fP^ K have unacceptably small p-values (p < 10~ 3 ). 

To better parametrize the /]| data, we are investigating alternative fit functions, including SU(2) 
S^PT. An improved treatment of the energy dependence of /jj would be especially desirable for 
adding data at higher momenta. (Cf. Sec. f|.) At the same time, f± dominates the desired form 
factor /+, and SU(3) S^PT fits for fP^ % and ff^ K are good. Any model dependence in f°-+* 
introduced by the #PT description of the energy dependence of fP^ n is probably small. We 
therefore proceed to compare the shape of the form factor fP_~^ n obtained from the fits shown in 
Fig. [I] to the shape measured by CLEO. 



3. Comparison of lattice QCD and experiment 

If lattice calculations are to be combined with experimental results, the lattice and experimen- 
tal results must be consistent. Experiments measure shapes of form factors but cannot fix their 
normalizations without the CKM matrix elements (Eq. (|1.1[)). For D — > K(n)£v, experiments and 
lattice calculations access the same q 2 region, and comparisons of the lattice and experimental form 
factor shapes provide stringent tests of lattice QCD. 

In Fig. |2| we overlay our calculated fP_~^ n , normalized to the point q 2 = 0.15 GeV 2 , with the 
same ratio from CLEO 0|. The orange (dark grey) error band is the statistical error obtained by 



including data from two coarse and three fine ensembles [13]. The yellow (light gray) error band 
is the statistical error obtained by including data from the 1 1 ensembles of Fig. |l| The curves are 
from SU(3) S#PT, and the errors are from 500 bootstraps. The errors scale as naively expected. 
The form factor shapes from CLEO and the lattice calculation agree well. 



4. Projected errors 

To conservatively estimate our errors, we begin with the error budget of our B — > 7t£v calcula- 
tion [Q]. At q 2 = naive scaling to the full QCD data set of Table [j] gives a statistical error of 4.2% 



5 



D — > K(%)tv form factors and \V cs u\\fwm lattice QCD 



Jon A. Bailey 



— 5 ensembles of full QCD lattice data (stat. errors only) 




0.5 I 1.5 2 2.5 3 

q (GeV 2 ) 

Figure 2: Overlay of the ratio 
f+^ K (q 2 )/f+^ K (q 2 ) from the lattice (curves 
and error bands) and CLEO (blue points) [S]. 
The larger lattice errors are from 5 ensembles, 
and the smaller are from 11 ensembles. The 
form factor shapes from CLEO and the lattice 
agree; the errors scale as expected. 




q (GeV 2 ) 

Figure 3: The form factors f+^ x {q 2 ) and 
fo^ K (l 2 ) an d their statistical errors. The curves 
are from S^PT, and the errors are from the Hes- 
sian. The percent errors are plotted against the 
right-hand axis. The larger errors at small q 2 prob- 
ably reflect excluded data at corresponding mo- 
menta on the finer, more chiral ensembles. 



and an error from the degenerate-mass vector-current normalization factors of 0.6%. Updating the 
heavy-quark and p -factor power counting estimates to account for the ultrafine data gives errors of 
2.5% and 0.5%, respectively. For the remaining errors we adopt our previous estimates [ 1C]. This 
leads to a total error of 6.1% for f^ K ( n > (fj). However, this estimate may be overly conservative. 
In Fig. H we plot the form factors (with the p-factors set to one) f®^ n (q 2 ) and f^^iq 2 ) and their 
statistical errors as functions of q 2 . The data are from the 11 ensembles of Fig. [I]. 

The errors at ^r^ ax are much smaller for fo than for / + because / + (/ ) is dominated by f± (/ji), 
and we have data for f x (f\\) for |p^| > 0.33 GeV o q 2 < 2.0 GeV 2 (|p^| >0^q 2 < q 2 ^). As q 2 
decreases, the errors reflect the addition of data and the hyperbolic behavior of the form factors; the 
errors in /+ (fo) grow approximately linearly in the region 1 .6 GeV 2 > q 2 > 0.8 GeV 2 (2.0 GeV 2 > 
q 2 ^ 0.6 GeV 2 ). The largest momentum of data in the fit is p^ = (1,1,0), which corresponds 
to q 2 G [0.43,1.58] GeV 2 on the ensembles with m l < 0.2m s , and to q 2 G [0.69,1.58] GeV 2 on 
the superfine ensembles. Without data points below q 2 = 0.43 GeV 2 on the finer, more chiral 
ensembles, the errors increase rapidly as q 2 decreases in the region 0.3 GeV 2 > q 2 > 0. 

Extrapolating the curve for the statistical error in / + to q 2 = 0, the error grows to about 3.6%, 
somewhat smaller than that expected from naive scaling to the entire full QCD data set. Including 
data at smaller q 2 would allow interpolation to q 2 = and might improve the errors significantly. 
We can appreciate the potential of the additional data by linearly extrapolating the curve for the 
error in / + for q 2 G [0.8, 1.6] GeV to q 2 = 0. The resulting expected error is about 2.0%. Adding 
a statistical error of 3.6% (2.0%) to our systematics yields a total error of 5.7% (4.8%). Corre- 
sponding error budgets are in Table f|| We conclude that including data at momenta greater than 
Pk(x) = (1,1,0) may improve the error in /+(0) to better than 5%. This prospect further motivates 
us to consider alternatives to SU(3) S^PT for describing the energy dependence of fu at small q 2 . 
(Cf. the last paragraph of Sec. |2[) 

Computations were carried out with resources provided by the US QCD Collaboration, Ar- 
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Stat. 


gn 


n 


m 




Kb 


Pn 


HQ 


Zy 


P 


FV 


Sys. 


Tot. 


(a) 


3.6 


2.9 


1.4 


0.3 


1.3 


0.2 


0.1 


2.5 


0.6 


0.5 


0.5 


4.4 


5.7 


(b) 


2.0 


2.9 


1.4 


0.3 


1.3 


0.2 


0.1 


2.5 


0.6 


0.5 


0.5 


4.4 


4.8 



Table 2: Projected error budgets for the form factors at q 2 = 0, assuming we (a) exclude data at momenta 
greater than = (1,1,0) and (b) include data at momenta greater than Pk(k) = (1, 1,0). Errors are due 
to limited statistics and the truncation of SjPT; uncertainties in the D*D% coupling, scale, average up-down 
quark mass, strange quark mass, and charm hopping parameter; momentum-dependent discretization effects 
of the light quarks and gluons; heavy-quark discretization effects; uncertainties in the matching factors Zy 
and p; and finite volume effects. The last two entries are the total systematics and total error. 
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